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Abstract 



We discuss numerical solutions of Einstein's field equation describing static, 
spherically symmetric conglomerations of a photon gas. These equations imply 
a back reaction of the metric on the energy density of the photon gas accord- 
ing to Tolman's equation. The 3-fold of solutions corresponds to a class of 
physically different solutions which is parameterized by only two quantities, 
e. g. mass and surface temperature. The energy density is typically concen- 
trated on a shell because the center contains a repelling singularity, which can, 
however, not be reached by timelike or null geodesies. The physical relevance 
] of these solutions is completely open, although their existence may raise some 

£f) • doubts w.r. to the stability of black holes. 

o 

On 

1 Introduction 

o 

The starting point of this investigation was the discussion of the Carnot-Bekenstein- 



^Jy process in the environment of a Schwarzschild black hole | HUS |. There it was 



assumed that the black hole is surrounded by a cloud of radiation with a local 



temperature according to Tolman's equation | TOL 



where is the usual Hawking temperature at r = oo. This equation also fol- 



lows for the equilibrium distribution of photons in relativistic kinetic theory [NEU] 



or from elementary thermodynamical gedanken experiments [ETUS]. These yield a 



modified formula for the efficiency of the Carnot-Bekenstein-process 



r, = i _ Wgoo(2) 
7W5oo(l) 

between two heat reservoirs at different height. In the equilibrium we have rj = 
which implies ([]]). All this holds as long as the back reaction of the radiation to the 
metric can be neglected. 

But, since goo( r ) = 1 — -y, at the Schwarzschild radius r = 1Z the temperature 
and the energy density diverges and cannot longer be regarded as a mere pertur- 
bation. Rather one would have to solve the field equation anew, this time allowing 
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for the energy-stress tensor of the photon gas to act as a source of gravitation. It is 
not at all clear whether the resulting metric will be a modified black hole in some 
sense. So we have the following problem which may be considered independently of 
the original motivation: 

Calculate the static, spherically symmetric metric of Einstein's field equations 
with the energy-stress tensor of a perfect fluid 

Tab = (p + P)u a u h + Pg ab (3) 

consisting of photons, i.e. 

p = 3P. (4) 

In section 2 the corresponding field equations are transformed into an autonomous 
two-dimensional system of differential equations which is discussed in section 3. In 
section 4 we investigate the metric of a photon star for r —* and study its global 
properties by means of numerical solutions. Physical characteristics like radius, mass 
and temperature are defined in section 5. Section 6 contains concluding remarks. 



2 Transformation of the field equations 



With except of eq. (|j) the problem stated above is just the well-known problem of 
constructing the interior solution of a star. We may choose coordinates (t, r, 9, 4>) 
such that the metric is given by 



ds 2 



-f(r)dt 2 + h(r)dr 2 + r 2 dtf 



(5) 



where / and h are unknown functions and dCl 2 is the surface element of a unit 
sphere. With respect to these coordinates the field equations boil down to a system 
of three coupled differential equations ( cf. [WAX] 6.2.3-6.2.5) 
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We recall that Tolman's equation [TOL] was derived under the same assumptions 
we made, except spherical symmetry. Hence we may adopt (||) or, equivalently, 



p(r) 



Pi 



f 2 (rV 



0) 



since 



P 



(10) 



from local statistical mechanics, where a is the Stefan-Boltzmann constant. We 
write 



Pi 



8vrG 



C 



(11) 
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and obtain from (||) and (0) 

C _ h! + h - 1 - 12 . 
/ 2 rh 2 hr 2 

3/ 2 rfh hr 2 ' 1 j 

Eq. dl) is identically satisfied, which confirms Tolman's result. Otherwise the 
system of differential equations would be overdetermined. 

The set of solutions 5 of the system (p^), ([i"3|) can be parameterized by 3 pa- 
rameters, e. g. C and two initial values for / and h. But S is invariant under the 
2-parameter group of scale transformations 

\2 

r h+ Ar, / ■-> ///, fcnli, C ^ -5C (14) 

A* 

Thus there exists only a 1-parameter family of solutions looking qualitatively differ- 
ent. 

An equivalent second order equation is obtained by eliminating h using 

h ~ C r 2 + 3f 2 (15) 

and inserting the derivative of (15) into (12): 

= 6r/ 2 + 6r 2 ff - 6/ 3 /' + 2r 3 f 2 
1 Cr s f + 3r/ 3 1 j 

For C = we have the well-known equations which lead to a 2-parameter family 
of Schwarzschild metrics. For C / we may scale every solution such that it 
becomes a solution with 

C = l. (17) 

The remaining subgroup of ( |i~4"D with A = fi may be used to simplify the differential 
equation (|l6|). We perform the transformation 

s = ln-, x = ^-, y = f(r)+x, (18) 
ro r 

where ro > is arbitrary. The resulting reduced equations read 

dx , , 

— = y-2x (19) 
as 

= ^-3x(x 2 -l)) 
(is x(3x 2 + 1) 

Note that this system is autonomous. The resulting symmetry s i— > s + so reflects 
the scale invariance of (]i~6| ) w.r. to the subgroup A = fi. Even if we cannot solve it 
exactly, it seems to be better accessible for intuition and developing approximation 
schemes. Any two different solutions of (^), (p0|) correspond to different similarity 
classes of S. 
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3 Discussion of the reduced equations 

We calculate some typical solutions of (|l9[), (pC|) numerically and display them as 
curves in the x — y-plane. 
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Figure 1: A selection of numerical solution curves of the 
reduced equation ([jl]), ( |20| ) together with the parabolic 
Schwarzschild approximation 

It is obvious that _ _ 

!BQ = yf, = (21) 

is a stable stationary point of (|l9|), (|20| ) which is an attractor of the whole open 
quadrant x > 0, y > 0. It follows by inverting the transformation ( |l8| ) that 

/(r) = y?r, %) = J (22) 

is an exact solution of (|l9|), ( f20| ) which is asymptotically approached by any other 
solution for r — > oo. It follows that the space-time of photon stars is not asymptot- 
ically flat. The other exact solution of (p^D, ( pC| ) 

y = 0, x = ae~ 2s (23) 

yields 

/(r) = -, h(r) = (24) 
r 

and is hence unphysical. A typical solution curve of (|l9|), ( ^0|) starts from x = +oo, 
y = 0, s = — oo and runs close to the solution ( |23"| ) until it reaches small values of x. 
Then according to the "x" in the denominator of (|20|), ^| increases rapidly and the 
curve is turned up towards the y-axis. Then it describes a parabolic-like bow and 
approaches the stationary point by a clockwise vortex. 

It is instructive to draw the general Schwarzschild solutions ( C = 0) 

fs(r) = a - - (25) 



4 



into the x — y-diagram. They are given by the family of parabolas 



x = y 

or 



(26) 



which approximate the solution curves of (19), (p0| ) having the same vertex at 



(27) 



In this way, for each solution (f(r),h(r)) of (|12|), (|13f ) we can define a unique 
SCHWARZSCHILD approximation (fs(r),h s (r)). 



4 Properties of the metric 

We now turn to the discussion of the solutions of (|l^ ) for f(r) which yield h(r) by 
©• 

To study the behaviour for smaller r we expand / into a Laurent series and insert 
this series into (|l6|). It turns out that the series starts with 

/(r) = -+£ + ■•■ (28) 
r 

in accordance with a Schwarzschild solution for C = 0. The coefficients A, 1? 
are left undetermined since they represent the two initial values for (|0]). The next 
terms are uniquely determined by (|l6l). It is straight forward to calculate the first, 
say, 20 terms by using a computer algebra software like MATHEMATICA. Here we 
only note down the first nonvanishing extra terms for / and h : 

A CB 

f(r) = - + B+ — + 0(r 5 ), (29) 
h~\r) = A + 1 __^_ r 4 + C?(r 5 ) . (30) 

For small r, / and h are also approximated by Schwarzschild solutions, but unlike 
the approximation discussed above, we have to choose A > in order to obtain a 
positive solution f(r) > 0. / cannot change its sign in a continuous way, since 
f( r o) = implies /"(ro) = oo by (|l6|). (According to the scale invariance (|i~4|), f(r) 
may be multiplied by —1, but this gives no physically different solution.) Thus we 
may state that for r — > the metric looks like that of a Schwarzschild black hole 
with negative mass, independent of C. For the geodesic motion close to r = we 
may thus adopt the effective potential of Schwarzschild theory ( |WAL| 6.3.15) 
(with M i — ^ — M): 

1 M I? Ml? , . 
V = K + K —+ + (31) 

2 r 2r r 



where 



(timelike geodesies) , . 

(null geodesies) 



It follows that r = can never be reached by particles or photons due to the 
infinite high potential barrier. Although curvature blows up for r — > 0, as in the 
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Schwarzschild case, the nature of the singularity is less harmful. We conjecture 
that it cannot be regarded as a "naked singularity" in whatever technical sense (see 
EAR] for details of the various definitions) and that Cauchy surfaces still exist in 



the spacetime given by (|l2|), (|l3|), if the line r = is excluded from the spacetime 
manifold. 

From a computational point of view, the singularity of f(r) at r = suggests to 



transform (16) into a differential equation for 

F(r) := rf(r) (33) 

and to re-transform to f(r) after a numerical solution for F has been obtained. We 
used the NDSolve-command of MATHEMATICA to produce the following numerical 
solutions: 




Figure 2: A typical numerical solution f(r),h(r) of the sys- 
tem (|l2j), (|l~3|). h(r) has its maximum at r = r$. The cor- 
responding Schwarzschild solution /s(r) with fs(ro) = 
is also displayed. 

A typical solution is shown in Fig. 2. Recall that for the Schwarzschild metric 
h{r) diverges at r = 1Z and f(1Z) = 0. For the solution of Fig. 2 h(r) has a 
relatively sharp maximum at tq and /(ro) is becoming small. For r > tq, f and h 
are comparable with their Schwarzschild approximations fs and hs, if r is not 
too large. For r < tq, f remains small within some shell r\ < r < tq and diverges for 
r -> according to @. By (g) this means that the energy density is concentrated 
within that shell and tq may be viewed as the "radius of the photon star". 

Other solutions with larger values of C show a more diffuse cloud of photons and 
a less sharp maximum of h, see Fig. 3 and 4. These solutions are all scaled to the 
same value of ro- 
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Figure 3: Numerical solutions h(r) for different C. The 
solutions are scaled such that they obtain their maximum 
at the same value tq = 1. 
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Figure 4: Numerical solutions f(r) for different C and the 
same scaling as in Fig. 3. 



5 Physical parameters of a photon star 

We have seen that the set of solutions S may be characterized by 3 parameters, 



e. g. A, B and C in p9|). From the analogy with the Schwarzschild case (C = 0) 
we expect that only a 2-parameter family represents physically different spacetimes. 
In the Schwarzschild case, one parameter is set to 1 by the choice of the units, and 
the remaining parameter 1Z distinguishes between black holes of different mass. More 
specifically, one postulates that the velocity of light, expressed by 4r approaches 1 
for r — > oo. The metric then obtains the form 

f S (r) = l-j, h s = f s \ (34) 
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In the case of the photon star, we cannot proceed in the same way, since the metric 
will not be asymptotically flat. But instead we may postulate the "gauge condition" 
that the Schwarzschild approximation of (/, h), defined in section 3, should obey 
condition (|34|). If this is not the case, one has to perform a suitable scale transfor- 
mation (0). In this way we obtain a two- fold of physically different solutions. 
We now consider physical parameters characterizing this two- fold of solutions. One 
could be the radius rg of the photon star defined above. 



In analogy to the Schwarzschild theory ( |WAL|1 6.2.7) we introduce the (gravita- 
tional) mass function 

2 

m(r):= r -(l-h- x {r))^. (35) 

In the domain where f(r) ~ fs(r) this is the "would-be- mass" of an equivalent 
black hole. In the domain r < tq the interpretation of ( |35| ) is not so obvious. As 
to be expected from the above discussion of the metric for r — > 0, it turns out that 
m(0) < 0. A typical mass function is shown in Fig. 5, where also the "proper mass" 

47r r i 

m p( r ) = — o~ / p(r)h(r)i dr (36) 
c Jo 

and the difference m p — m is displayed. 




Figure 5: A typical numerical solution of the mass functions 
m(r),m p (r) and the difference m p (r) — m{r). 

It may be, as in this case, that the majority of the photons are "hidden" by the 
apparent negative mass in the center with respect to gravitation. 
Nevertheless, we could use 

m Q = m(r ) (37) 

as a further physical parameter characterizing a photon star. Generally, by (^) and 
K r o) > f = Hm?--^ h(r) 

^M s (r ) <m < M s (r ), (38) 
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where 



M s {r ) 



2G 



(39) 



If h(r ) » 1 we have m « M s (r ), as in the Schwarzschild case. 
As another physical parameter we consider the "surface temperature" 

>(r )^ 1/4 



'-o := 



T(r ) 



(40) 



If we have only a two-fold of physically different solutions, as we claimed above, To 
should be a function of r$ and rriQ. Indeed, (|l2| ) together with h'(ro) = shows that 
f(ro) is a function of tq and h(ro), hence of ro and itiq. Then, by (40) and @, also 
To depends only on ro and tuq. Since a depends on h, the result can be conveniently 
expressed by using Planck units, indicated by a subscript P: 



To 



T P (2tt) 3 / 4 



m ( L P Y 
M P \r J 



1/4 



(41) 



To give a numerical example, if we take as the mass of the sun, ro as the 
corresponding Schwarzschild radius, we obtain T « 4 • 10 12 K. This would cor- 
respond to a very hard gamma radiation with a wavelength A ~ 10 -15 m. The 
Hawking temperature of this example is T# ~ 10 _8 K, since T# ~ m _1 whereas 



'o 



m, 



-1/2 



for m ~ Ms(tq). 



6 Conclusion 

It is difficult to assess the physical relevance of our findings. But one point seems to 
be clear: The global character of the solutions (/, h) is completely different from the 
Schwarzschild approximations (fs,hs), no matter how small C is. So the class 
S of solutions of (|I2"D, ( |l"3| ) does not depend continuously on C in any reasonable 
sense. Any small amount of radiation will destroy the event horizon, at least if an 
equilibrium is approached. If this would also be the case in simulations of the birth 
of "black holes" by collapsing matter, then they would never be born, and perhaps 
won't exist at all. 
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